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ABSTRACT 


An analytical investigation was carried out to determine the 
design conditions, for a ducted propeller unit operating in water, 
providing optimum suppression of cavitation when the forward speed, 
depth and body drag are given. The selection charts were constructed 
in detail for propellers using the NACA 65 series blades. The results 
indicate the pumps of generally higher flow rates are more desirable 


than those that have deen employed in the few units actually constructed. 
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LIST OF SYMBOLS 


cross-sectional area at entrance to guide vanes (£7) 
drag coefficient of propelled vehicle 

drag (ib) 

thrust (1b) 


cavitation number associated with minimum operating depth, 


at the inception of cavitation 

atmospheric pressure (inane a) 

local static pressure (ibe / ft”) 

vapor pressure of water (ibs /ft”) 

total pressure (lbs / ft”) 

radius of propeller from axis of rotation 

pressure coefficient across fan blades 

propelled vehicle cross-sectional area ct”) 

velocity relative to rotating blade (ft/sec) 

chord length 

minimum operating depth (ft salt water) 

tangential blade spacing (ft) 

axial velocity, absolute (ft/sec) 

tangential velocity, absolute (ft/sec) 

number of blades 

angle of attack between flow direction and fan blade chord 
flow angle between flow direction and cascade axis 

fan blade stagger angle between cascade axis and blade chord 
propulsive efficiency 


turning angle, the difference between flow angle at entrance 


to fan blade and exit from fan blade 





Na 


drag-area coefficient, propelled vehicle cross-sectional area 


divided by duct cross-sectional area, times the drag coefficient 


diffusion coefficient, absolute axial velocity at inlet to guide 


vanes divided by free stream velocity 
density of sea water (lbs /ft?) 
solidity, chord of blades divided by tangential spacing 


flow coefficient, rotating velocity divided by free stream 


velocity 


turning coefficient, tangential turning velocity divided by free 


stream velocity 


J angular rotating velocity 
subscripts 

O free stream 

upstream of guide vanes 

2 upstream of rotating fan blades 

3 downstream of rotating fan blades 


downstream of nozzle 








I. INTRODUCTION 

The ducted propeller is of interest in hydrodynamic applications 
largely because it offers the possibilities of either high loading at low 
speeds or cavitation suppression at higher speedz. The present work 
igs concerned primarily with utility of the ducted propeller to suppress 
cavitation on the propeller blades as a means of achieving reasonably 
efficient and quiet operation at high speeds. The several actual attempts 
made to achieve this end have not resulted in inarked success largely 
because the principle was applied to problems where the speed, depth 
or propeller design were inappropriate. This difficulty is under- 
standable because, at moderate torward speeds, the selection of 
propeller blading, rotative speed and diffuser design must be made in 
a rather particular manner to achieve the potential advantage of the 
ducted propeller. Furthermore there existed no method for selection 
of the appropriate configuration to satisfy given requirements. It is the 
purpose of the present analysis to develop a method for finding optimum 
configurations for suppression of blade cavitation and to carry out 
design charts in detail for a particular class of propeller blading. 

The conceptual model employed is an axial flow ducted fan with 
cuide vanes ahead of the rotating fan blades. Possible physical 
arrangements are shown in figure la, bandc. The flow entering the 
hydrojet of figure ia and lb is that flow from the wake of the vehicle 
being propelled. The flow entering the hydrojet of figure le is the 
free stream flow. This investigation was limited to considerations of 
the free stream flow entering the hydrojet. Cavitation was assumed 


to occur on the fan blade; no consideration was piven to the possible 
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occurrence of cavitation due to the geometry of the diffuser and duct or 
of the body being propelled. Although it is known that cavitation may 
occur first at the blade tip or hub, high noise output, loss of efficiency 
and structural damage are largely associated with cavitation of the 
blade surface. 

In the development of the analysis, the free stream flow velocity 
U,. was considerea to >e ideally diffused to a lower velocity u, at the 
entrance to the guide vanes. The flow through the cuide vanes and 
rotating propeller blades was analywed using airfoil cascade theory. 
The flow was treated at a typical radius R from the rotating axis and 
is thus approximately applicable to any section of the propeller blade. 
The guide vances were considered to impart a tangential velocity, nw , 
to the flow entering the rotaiing blades, tne rotating Llades were 
considered to turn the flow back to an axial throughflow, and to 
increase the total pressure by AP,. IK was assumed that the total 
pressure rise was constant from hub to tip, that the axial flow velocity 
was constant, and that the radial velocity was negligible. The flow 
downstream of the retatin; fan blades was then discharged through a 
nozzle to an exit velocity uy and an exit pressure py , cgual to the 
free stream pressure p,. Cavitation occurred on the suction side 
of the rotating blade section when the local static pressure p, dropped 
to a certain predetermined value p, . 

The thrust of the unit was determined by elemcntary application 
of the momentur. theorem to the flow upstream and downstream of the 


ducted propaller. Since steady motion was assumed, the thrust must 
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equal the body drag and hence propeller requirements were related to 
the drag characteristics of the body. 

Performance data for the rotor blades was extracted from 
references (1) and (2) for the NACA 65 series airfoil incascade. These 
results were employed to obtain a typical set of selection charts, 
anslicable only to these blades, to illustrate use of the analysis to select 
appropriate blade sections for a given forward speed, depth of operation 


and body to be pronelled. 





ae 
Il. DEVELOPMENT OF THE ANALYSIS 


The thrust of a ducted propeller is derived frorn the change in 
momentum of the fluid flow thru the unit. This change in momentum 
ig accomplished by the transfer of energy between the fluid and the 
rotating propeller blades resulting in a total pressure increase which 
is converted to kinetic energy in the nozzle. The mass flow thru the 
duct is described by OU A,, where ¢ is the density of sea water, U, 
is the absolute axial flow velocity upstream of the guide vanes, and A, 
is the cross-sectional area of fluid flow at the inlet to the guide vanes. 
By continuity for fluid flow and since ? is constant for an incom- 
pressible fluid, the mass flow upstream and downstream of the hydro- 
jet is also GU A,. The change in momentum thru the ducted pro- 
peller is therefore the mass flow times the change in velocity between 
upstream and downstream states, and is equal to the thrust of the hydro- 


jet 
—_— © us Ai (uy Sa (Z 1) 


This thrust must equal the drag of the unit being propelled in 


order to maintain a forward speed of Uo. Drag is expressed as 
/ 
D=+ euGS (2.2) 


where Co is the drag coefficient defined so as to include, in addition 
to usual aerodynamic drag, the pressure forces over the streamlines 
separating the flow passing through the pump from that passing around 
it. S" is the cross-sectional area of the body being propelled. From 


equations 2.1 and 2.2 we have 











(2) + (2. 3) 


Propulsion efficiency is defined as the thrust work divided by 


the work input and is expressed as 


gtr Ar (Ys ~ Uo) Wo (2. 4) 


Ne ss ie = Teac Ay 
zZ 


which reduces to 


t) 


Ne Tf uy ico 
a ( Uo + | 
From equations 2.3 and 2.5 it is seen that two major dimen- 
sionless parameters appear that are important to efficiency and to 
the size requirements of the hydrojet relative to the size of the vehicle. 


These are defined as 


Zr C e ne 
RT A ay, ia) 





where ) is a combined drag and cross-sectional area coefficient 
and 7T\a is a diffusion coefficient from free stream to the guide vanes. 
The interrelations of A , 7\q and Np are shown in figure 3. 

To analyze the guide vanes and rotating fan, the airfoil cascade 
theory was applied. An imaginary cylindrical cross-section of the 
blades is taken with axis coinciding with the axis of rotation. The 


imaginary cylinder is assumed to have radial extent dR anda radius R 





and then is imagined to be unrolled. A row of parallel blade cross- 

sectional cuts results which extends on to infinity if the pattern is re- 

peated indefinitely. The flow problem then is reduced from a three- 

dimensional to a two-dimensional problem. The rotational velocity of 

the rotor blades is represented by a tangential velocity wR perpendic- 

ular to the cascade axis. An observer can assume that he is fixed in 

a coordinate system relative to the rotating blade sections or in an 

absolute system with respect to the guide vanes. The difference in 

the two coordinate system then is only in the rotating velocity, WR. 
The relative spacing of blade sections is described by the 


solidity, O , and 


oS c/t (2.7) 


where C is the blade section chord length and t is the tangential 
Spacing equivalent to 27(R_ divided by the numoer of blades. The 
solidity parameter is one of the important design quantities which 
effects the performance of the airfoil sections and the occurrence of 
Cavitation. 

In figure 2 is shown a typical row of guide vanes and rotor blade 
sections obtained by the unrolling of a cylindrical cross-section. 
Reference stations are designated as (0) to indicated tree stream flow, 
(1) to indicate the flow entering the guide vanes, (2) as exit frorn the 
guide vanes and entrance to the rotating fan, (3) as exit from the 
rotating fan, and (4) as the exit from the nozzle. The angle between 
the blade section chord and cascade axis is designated the stagger 


angle, JS . The direction of flow relative to the rotating blade sections 
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is established by the inlet flow angle, (2, measured from the cascade 
axis to the direction of flow. Similarly, the exit flow angle (4; , is 
measured from the cascade axis to the direction of flow relative to and 
leaving the rotating blades. The difference between the inlet flow angle 
and the exit flow angle is the turning angle, 6 , thru which the flow is 
turned by the rotating blades. 

The flow through the ducted propeller, absolute and relative to 
the rotating blades, can be represented by a velocity triangle as shown 
in figure 2. The axial flow velocity was assumed constant therefore 
U, = Uy . The flow velocity relative to the rotating fans is obtained 
by suotracting the rotating velocity, wR , from the absolute velocity, 
U, , entering the blades sections, giving the inlet flow velocity, We P 
relative to the rotating blades. The relative flow velocity leaving the 
fanis W3 , and is converted to an absolute velocity, u, , by adding 
the rotating velocity, wR, to the relative velocity. The tangential 
velocity, Ar , imparted to the fluid by the rotating blades was assumed 
to be equal to the tangential velocity imparted by the guide vanes, 
resulting in the simplified velocity triangles as given in figure 2. 

To describe the pressure distribution across the rotating blade 
section, the pressure coefficient, S , is used as given in references 


(1) and (2) where 


& Par — Px 
$ i ewe (2. &) 


the total pressure entering the rotating blade section, minus the local 
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static pressure at any point across the blade section, divided by the 
inlet dynarnic pressure. When the local pressure, R , is equal to 
some predetermined pressure FP, , cavitation is assurned to occur. 
The value of F, may be taken as the local vapor pressure of the 
water or a more suitable figure if knowledge of local conditions 
permits it. Thus for given inlet conditions there will be a maximum 
value of the pressure coefficient, 5 nax , fOr Operation without 
cavitation. To relate the free stream conditions to Snax , a 
cavitation number \ , is employed where 


K = foo (2.9) 


Zz Pe! 

the free stream static pressure minus the cavitation pressure of sea 
water, divided by the free stream dynamic pressure. It will be shown 
that the value of K determines the minimum Operating depth d 
corresponding to P, for a given operating velocity, U, . Part of the 
problern of this analysis is to relate the cavitation number K , to the 
pressure coefficient, Spay , where Smax is known frorn experi- 
mental data. 

Two other dimensionless parameters employed for this analysis 


are the flow coefficient Q , and turning coefficient WY where 


QP = wR (2. 10) 


Uo 


ys + (2.11) 








thus relating the rotational velocity and the turning velocity to the free 
stream conditions. 

By the assumption that no losses occur, the total pressure rise, 
AP, across the propeller can be obtained by subtracting the free stream 


total pressure from the exit total pressure so that 


OG Pig dala (2.602) 


The total pressure is equal to the sum of the static and dynamic 


pressures so that 


Par = Pao + 5 7 Ud (2203) 
and 
z - 
fomeaviito. aes g Uo (2. 14) 
Since the static pressures at entrance and exit are assumed to 
be equal, 
AP, =t oe Se) (2. 15) 
1 $ wes : 


However, since no losses occur and the axial flow velocity was 


assumed constant, we obtain 


fete ry eee oe (2. 16) 


Pe Goes ek (255i) 








and 


AP, = P, —P, (2213) 


Thus, Bernoulli's equation relative to the rotating Slades can be applied 
to obtain the total pressure rise across the hydroject. Therefore, 


relative to the rotating blades 


é 
Pot aeW, = At 4 pws (2. 19) 


From the velocity triangle of figure 2, it is seen that the relative 


velocities are 


ry 


oe 2 ee (2. 20) 


si 


i SSS (eR (2.21) 


where again the axial velocities U, and UU, are equal. Relative to 
the guide vane, Bernoulli's equation would be 


7 


Pttous = Pi +h eU2 (2. 22) 


where 


rs 2 t 
a > omer (22a) 


By combining equations 2.18 thru 2.23 Euler's turbine equation is 


ootained as follows: 


DP SP, = Pe ee (2. 24) 
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Therefore from equations 2.15 and 2. 24 we obtain 


ca 


bh = ewRoe= teu [[) -jj (2. 25) 


& 


or in terms of the flow coefficient Q and turning coeificient a 


roe 
CGY = (Ss) — | (2. 26) 


When thrust and drag were equated it was found, equation 2.3, that 


Md gD +1 (2. 27) 
Therefore 
ower AL + Hr elAl a) (2. 28) 


This important relationship matches the operating characteristics of 

the ducted propeller, D and U , with the drag of the vehicle being 

propelled, the relative size of the vehicle and duct cross-section, and 

the performance of the entrance diffusor. Ina sense, the product 

is a measure of the power requirements of the ducted propeller. 
Associated with the rotational speed and turning velocity is 

the turning angle, © , which is the difference in the flow directions at 

the entrance and exit relative to the rotating blades, therefore 


tayo = tan Ee — tan 83 2.2 
[4 Cai oreean, 3, ( m 


But from the velocity triangle of figure 2 it is seen that 





Jpn. 


2 6h vik 
Se as (2. 36) 
tan Gy = rs (2. 31) 


By substitution and by dividing numerator and denominator of equation 
2.29 by the square of the free stream velocity, the turning angle can 


be determined to be 


Q@= Tle W 
tan Dig+yvi+ Te (2562) 


From: equation 2.32 a value of © can be calculated which then must 

‘e metched with experimental performance data for the particular 

type airfoils being employed for the propeller blades. Associated with 
this turning angle will be a pressure coefficient, S , which can be 
allowed to have a maximum value when the local pressure across the 
blades is equal to the pressure Pp, . 


From the definition of the pressure coefficient we have that 


oe = 
G2 ee eee (2233) 
+ oW 4 OW 
z 5 “ Z y a 
But 
+ 2 2 
R=Rttelue-—(wra) ] (2. 34) 
ee See ; > <i L 
therefore by dividing numerator and cenoiminator vy FP Us 


we obtain 
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R, _ Py 
Fue! 





+ [i-T -¥] 
[pt] + The 


7a Nag 


a= (2. 35) 


When the local pressure equals the pressure Py we have Snax and since 
the cavitation number M was defined as 


‘e ee 


oe 


ie ys 
x § He 
we can determine the cavitation number as a function of the ducted 


propeller performance and operating characteristics 


Ke (Sue NOW) +MY +EE MN FV -1 (2. 36) 


From the definition of the cavitation number it is seen that for a given 
operating velocity U. the minimum operating depth Gd is determined 
since the free stream static pressure ee is cquivalent to the pressure 


at the depth d plus the atmospheric pressure. The expression relating 


K » Uo and d , assuming the vapor pressure of water is zero, 15 
de ue — & (2. 37) 
“sg Fg 


and is presented as a nomograph in figure 19. Therefore, frorn equations 
2.36 and 2.37 the operating depth and velocity can be determined as a 
function of the performance of the rotating blades, the diffusion and 
power requirements for the ducted propeller. To operate at shallow 
depths with high velocity it is seen that the pressure coefficient should 


have a minimum value for the maximum suppression of cavitation. 
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This in turn requires that the pressure coefficient should have as 
small a value possible which is equivalent to saying that the velocities 
across the blade sections must be minimized. Therefore, to suppress 
cavitation there will be optimum values for all of the dimensionless 
parameters discussed which will match with the actual performance of 
the blade sections as given in references (1) and (2) for airfoils in 
cascade. 

In this section, the important relationships have been derived. 
In the next section the computational procedure carried out to solve 
the relationships and to match the solutions with the experimental 
performance data for the 65 series airfoil sections in cascade will be 


given. A summary of the important relations are as follows: 


= oe IS 
oy A Late 26 | (2. 28) 
| 
a en ee (2.5) 
ama * | 


‘e on = eel Li eee 2182 
a O[O+y] + hy ; 


K = (Sux [(O+y) +A] + [Ag+ v*- 1] (2. 36) 














Lag 
ll. COMPUTATIONAL PROCEDURE 


From equations 2.5 and 2. 28 it is seen that for fixed diffusion 
coefficient Ne and drag-area coefficient ) , the propulsive efficiency 
Np, and the product of the flow coefficient Q and turning coefficient 
are fixed. Then using equation 2.32, the turning angle © can be 
determined for a range of values of ( for a series of constant values 
of [\4y and dh . This was done for values of 7\4 between .4 and 
1.1 and for values of A between .25 and .75. It was necessary then 
to match the computed values of the turning angle with the performance 
data for 65 series airfoils in cascade in order to determine the 
pressure coefficient, 5y,.. 

The pressure coefficient S yj, is a function of the solidity, Oo , 
angle of attack, ~ , stagger angle, 5 , camber and thickness dis- 
tribution. The data presented in reference (i) gives S and @ asa 


- 


function of o& for values of solidity equal to .5, .75, 1.00, 1.25 
and 1.50, and for values of flow inlet angle, G,, equal to 30°, 45°, 
60° and 70°. The data is also presented for various airfoils of the 
65 series with same thickness distributions but different values of 
camber. In reference (2), data is presented for airfoil sections with 
the camber reversed in an attempt to reduce velocities across the 
blade sections. 

Fixing the inlet flow angle (. , and varying the angle of attack 
ox » requires that the stagger angle be varied. However for a fixed 


blade row the stagger angle must remain fixed requiring that the inlet 


flow angle vary with the angle of attack. Therefore to use the data 
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presented in references (}) and (2) it was necessary to cross plot the 
data to obtain data for fixed airfoil geometry. The cross plotting was 
carried out for stagger angles of 45° and 55°, Figure 5 isan 
example of the results for the 65 4:0 airfoil in cascade with a solidity 
of ..0, stagger angle 55°, and with Sugand © plotted vs « 

The computed values of turning angle, S& , were then entered 


into the cross plots to determine the maximum value of the pressure 


coefficient, Snax , corresponding to the values of Ny . A . @ 
and W . Having Snax . the cavitation number K was calculated and 
plotted. Figure 10 shows an example of sucha plot for A - .5, 


s. = 55°. Oo - 1.0 where the curves represent values of K and Yy for 
fixed values of Tig. A cross plot was then made for constant values 
of the cavitation number MK as a function of T\y and ). Figure 13 1S 
an example of the contour plots of K obtained and since a minimum 
value of K is required for the maximum suppression of cavitation 
the design point for the given conditions is the center of the contour 


representing the minimum value of the cavitation number. 
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IV. DISCUSSION OF ANALYTICAL RESULTS 
To determine the effect of the solidity, three cavitation number 


contour plots were made for values of o = 1.0, 1.25 and 1.50 where 


/ 
Co 


the stagger angle was fixed at 55° and A= was fixed at .5. 
Figures 13, 15 and 16 are the result and clearly show that solidity 
should have a value of approximately 1.25. Similar plots for values of 
O- equal to .75 and .5 were not made because by inspection of the Spa 
and 8 vs & plots (not shown) the values of Smax are somewhat 
higher than for the values shown in figures 5, 6 and 7. 

The effect of stagger angle, which is another important design 
parameter, can be determined by inspection of 5,,and @ plots. The 
values of Sux for stagger angles less than 45° are larger near design 
point (minimum Snax ) but will give a larger spread of values of Pp 
for minimum values of K . For stagger angles greater than 55° the 
values of Smax are again larger and also the spread of values of Y 
igs more We In figures il and 13 a comparison for stagger angles 
of 45° and 55° can be made with all other conditions remaining fixed. 
The increase to stagger angle, S) ,» from 45° to 55° has shifted the 
design point slightly to the right, decreased the minimum value of the 
cavitation number, and increased the area inclosed by values of 
constant Ki . Any further increase to the stagger angle can be 
expected to shift the design further to the right and increase the minimum 
values for ( . The physical significance of shifting the design point to 
the right is that larger values of { can be employed which means that 


the rotational velocity can be increased for a given value of U, , and 
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that the loading of the blades is decreased by a decrease in the angle 
of attack. The optimum value for the stagger angle appears to be 
approximately 55° for the airfoils considered. 

The effect of drag-area coefficient, A , 1S quite pronounced 
in shifting the design point. Calculations were made and plotted for 
values of A =.25, .5G and.75. The value of A = 0.25 is generally 
appropriate for a good torpedo shape. The resulting cavitation number 
contour plots are shown in figures 12, 13 and 14. In addition, constant 
values of K are cross plotted to figure 3 to show the effect of \ on K 
as well as propulsive efficiency, Vp . From the figures it is seen 
that increasing 4 shifts the design point to the right, increasing Pp , 
and up, increasing the design value for the diffusion coefficient Ta. 
It also radically increases the minimum values of N_ that can be 
obtained, thus restricting the operating velocity and minimum depth of 
Operation. Therefore, for the maximum suppression of cavitation, 
the drag-area coefficient, A. , should be as small as possible. This 
requires a small value at the drag coefficient, Ceo , anda small value 
for the ratio of cross-sectional area of unit being propelled to the cross- 
section area of the duct. If the drag at the unit is high, a comparatively 
large ducted propeller will be required. The value of therefore 
has a marked effect in determining the maximum size of a vehicle for 
propulsion by ducted propeller. 

To reduce the problem of obtaining diffusion ratios other than 
1,0, a maximum value of A is desirable since the design value of Ta 


increases with A . In order to obtain a larger range of values for 
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the flow coefficient QM figures 12 through 14 indicate that a minimum 
value for A is required. Since % ~& for given speed and angular 
velocity, large range in YP is desirable to permit a large hub ratio 
and hence a small outer diameter. From figure 3 it is evident that 
decreasing » will increase the propulsive efficiency Np. In view 
of the above considerations, alow value for A is again more 
desirable than a high one. 

The effect of increase in camber of the blade sections is to 
increase the velocities acrogs the blade surfaces and thus increase the 
pressure coefficient, Sm. Calculations and plots were made for the 
65-410 and 65-810 airfoil sections for comparison where stagger angle 
and A were fixed. The results are shown in figures 13 and 17 and 
show that increase in camber from 4 °/o to 8 °/o shifts the design to 
the left, decreasing the design value for Y . The minimum value for 
the cavitation number, K , has not been appreciably affected but the 
areas enclosed by constant values of K have been reduced considerably, 
reducing the range of Y values available for a given design value for 


K . The vlots of and @ vs oc for airfoil sections with greater 


Sn 
cambers show that Snaxincreases considerably and that the range of 
useful values of Q will be reduced further. 

In figure 18, the resulting contour plot for the 65- (8A, Ia,)IO 
airfoil sections is shown to point out the effect of reversing the camber 
distribution to reduce local velocity ratios. Although this decreases 


the minimum value of K , the range of values for Q at fixed 7Ta , 


is decreased severely. Thus it does not appear advantageous to reduce 








r 


ldcal velocities over blade suriice by this technique. 
To illustrate how data given in the figures can be employee to 


2. 


determine optimum design conditions for the maxi:nurs suppression of 
cavitation, assume thet a unit is to ve propelled at 2 ,elocity oi 73 ft/sec 
(25 <nots) and is to operate at a depth of ter feet. Cron: the previous 
discussion it was seen that a minimur: value of the drug-area coefficient 
is required. Therefore, design consideratiors for the unit will .e for 
as low a drag coefficient, Cp , as possible, and for as favoraile a 
crosa-sectional area ratio as possible consistant with size, weight and 
space availavle to the unit being propelled. Assume tnat the best value 
of AX that can be obtained iy 6.5. From fipure i9 it is found that the 
meximum allowable value for the cavitation nu.abecr is approximately 
0.5. £0r any value greater than ¢.5, cavitation can de expected to 
occur. For A = 0.5, figure 15 for the 65-116 clade 
to have tne largest area enclosed by a cavitation number of §.5. It 
gives a large range for the values of QM and requires a Cciffusion 
veificient, T\g , of 9.61. The solidity and stayyer angle in this case 
are 1.25 and 55° respectively. The center of the contour representing 
the mini:num values of K would we designated as the design »oint 
for the mean radius R between hub and tip. This value for Y = 1.15 
then would fix the optimum angular velocity w and thus fixes the 
values for QP for changes in radius RK from: hul: to tip. The ratio 
of the maximum and mininiurn flow coefficient contained within the 
boundaries for K = 6.5 and T\y = 0.61 is approximately 1.3. This 


indicates that a relatively large size nus will De necessary, the actual 





Dm 


size being fixed by the value of cross-sectional area decided upon in 
fixing A. 

From figure 3, the propulsion efficiency can be determined. 
Then from the propulsive efficiency and drag coefficient, the power 
requirements can be determined by replacing the thrust work by drag 
and solving for work input in equation 2.4. Thus a design criterion 
has been established for the maximum suppression of cavitation. 

If optimum diffusion coefficient values could not be obtained, if 
greater velocities were required, or if the value of the drag-area 
coefficient A were unfavorable, it would be necessary to increase the 


minimum depth of operation and hence raise the free stream pressure 


_ 











area ic 


V. CONCLUDING REMARKS 

An attempt has been made to determine what design and operating 
criteria are necessary to select a ducted propeller propulsion unit with 
the maxirnum suppression of cavitation. The detailed selection charts 
were limited to the NACA 65-series airfoils for propelier blades. The 
PPeiqus, however, is general. 

The most sensitive parameter appears to be the drag-area 
coefficient, . This quantity should have as low a value as possible 
in order to increase the range of propeller operation without cavitation. 
The diffusion coefficient J\yis also a significant parameter and it is 
an interesting result that an optimum retardation of the stream exists; 
the stream may be slowed down too much as well as too little. 

The flow coefficient parameter, g , was found to have a rather 
low optimum value, requiring relatively slow rotating velocities or 
through flow velocities much larger than conventional. In addition it 
was found that the range of values for g without cavitation will be some- 
what limited. The design criteria for {) are largely determined by 
the blade sections being employed for the propeller, requiring low 
velocities across the blades with high loading of the blades. 

For the blade sections investigated, it is found that there are 
optimum values for the solidity O , the stagger angle, s , and the 
camber. The optimum values are those that will provide minimum 
values of the maximum pressure coefficient, one , over as largea 
range of angles of attack as possible, and with limit on the maximum 


values of the turning angle. 
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Figure 1, 


EA 


a) Wake flow entering hydrojet. 


b) Wake flow entering hydrojet driven by conventional 
propeller shaft. 





c) Free stream flow entering hydrojet 


Schematic of poss 


ible hydrojet propulsion configurations. 
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Figure 2. Schematic of guide vanes and hydrojet blade 
sections in cascade and the associated velocity 


triangle. 
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Figure 3. Diffusion coefficient, /\g , versus propulsion 
efficiency np, for constant drag-area coefficient, A . 
Constant values of K are shown for the combination 
o =1.25 , 5 =55° and NACA 65-410 blade section. 
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Figure 11. Diffusion coefficient J\4 vs. flow coefficient Q for 
constant values of cavitation number K for the com- 
bination A=.5, O0=1.00, J = 45° and NACA 
65-410 blade section. 
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Diffusion coefficient Tla vs flow coefficient Q for 
constant values of cavitation number K for the com- 
bination A =.25, O = 1.0, f =55° and-NACA ~ 
65-410 blade section. 
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Figure 13. Diffusion coefficient 7\4 vs flow coefficient PY for 





constant values of cavitation number K for the 
combination \ =.5, © = 1.00, JS = 55° and NACA 
65-410 blade section. 
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Figure 14. Diffusion coefficient T\a vs flow coefficient QP at 
constant values of cavitation number KM for the 
combination A = .75, o = 1.00, § = 55° and NACA 
65-410 blade section. 





Figure 15. 


-39- 





Diffusion coefficient Tlavs flow coefficient QP for 
constant values of cavitation number \ for the com- 
bination A =.5, oO =1.25, J = 55° and NACA 
65-410 blade section. 





-40- 


NCSL 
| eM 
monn) 


.6 ars. Lao) Q Le 1.4 L276 Lae 





.@ °) 


n 


fas 


Figure 16. Diffusion coefficient Tia vs flow coefficient YP for 
constant values of cavitation number K for the 
combination \ =.5, 0 =1.5, 4 =55° and NACA 
65-410 blade section. 
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Figure 17. Diffusion coefficient T\avs flow coefficient Y for constant 


6 8 Toko D 


values of cavitation number K for the combination A =. 5 


o =1,.00, J =55° and NACA blade section 65-810. 
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Figure 18. Diffusion coefficient T\a vs flow coefficient Y for 


constant values of cavitation number Y“ for the com- 
bination \ =.5, ¢ =1.00, S =55° and NACA 
65- (8A> Igyh)10 blade section. 
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Figure 19. Nomograph relating the cavitation number K , free 


stream velocity U., and depth of operation d. 



































